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Abstract Escherichia coli is a motile bacterium that moves up a chemoattrac- 
tant gradient by performing a biased random walk composed of alternating runs 
and tumbles. This paper presents calculations of the chemotactic drift velocity Vd 
(the mean velocity up the chemoattractant gradient) of an E. coli cell performing 
chemotaxis in a uniform, steady shear flow, with a weak chemoattractant gradient 
at right angles to the flow. Extending earlier models, a combined analytic and nu- 
merical approach is used to assess the effect of several complications, namely (i) 
a cell cannot detect a chemoattractant gradient directly but rather makes temporal 
comparisons of chemoattractant concentration, (ii) the tumbles exhibit persistence 
of direction, meaning that the swimming directions before and after a tumble are 
correlated, (iii) the cell suffers random re-orientations due to rotational Brownian 
motion, and (iv) the non-spherical shape of the cell affects the way that it is rotated 
by the shear flow. These complications influence the dependence of on the shear 
rate y. When they are all included, it is found that (a) shear disrupts chemotaxis 
and shear rates beyond 2s"' cause the cell to swim down the chemoattractant 
gradient rather than up it, (b) in terms of maximising drift velocity, persistence of 
direction is advantageous in a quiescent fluid but disadvantageous in a shear flow, 
and (c) a more elongated body shape is advantageous in performing chemotaxis 
in a strong shear flow. 
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1 Introduction 

Free swimming bacteria are present in many naturally occurring aqueous environ- 
ments including animal intestines (enteric bacteria) and the open ocean (marine 
bacteria). In the ocean, motility could influence bacterial ecology and the role 
of bacteria in oceanic biogeochemistry [19]. Thus there is an oceanographic mo- 
tivation to understand how motile bacteria behave in a weakly turbulent, sheared 
environment. Bacteria are also known to adhere to each other and/or other surfaces 
and form biofilms 1 13]. To understand the formation of biofilms in the presence 
of a fluid flow it may be helpful to model the motion of free swimming bacteria in 
the sheared region adjacent to the biofilm. 

In this paper we focus on modelling Escherichia coli, which is the most well 
studied motile bacterium. As described by Berg |@], E. coli is a common enteric 
bacterium with a rod shaped body w 1 /im in diameter and w 2 jUm long. A typical 
E. coli cell has w 6 left-handed helical flagella emerging from random points on 
the sides of its body, each of them extending several body-lengths into the sur- 
rounding fluid. The flagella are powered by reversible rotary motors. When all the 
flagella spin counter-clockwise, they form a synchronous bundle and propel the 
cell forward in an almost-straight 'run'. When one or more flagella spin clock- 
wise, the flagellar bundle comes apart and the cell swims in a highly erratic 'tum- 
ble', with little net displacement but a large change in direction. Tumbles exhibit 
'persistence of direction', meaning that the swimming directions before and after 
a tumble are correlated, with a mean angle between the new and previous swim- 
ming directions of w 62° |8]. Hereafter, for brevity, 'persistence of direction' shall 
be referred to simply as 'persistence' . 

Even during a run, the cell exhibits small, random changes in direction, and 
these have been attributed to rotational Brownian motion {i.e. thermal collisions 
with molecules in the surrounding fluid). Berg [6] estimated the coefficient of 
Brownian rotation to be Z)r w 0.062 radians^s"' for an E. coli cell swimming 
in a fluid of viscosity 2.7 cp at 32°C, and this is consistent with experimental 
observations |8]. However, this estimate was based on treating the E. coli cell as 
a sphere of diameter 2 /Jm, and when one takes into account the stabilising effect 
of the flagellar bundle, the theoretically predicted coefficient of Brownian rotation 
is an order of magnitude smaller than the observed rotational diffusivity [17], so 
it seems likely that the observed diffusivity is in fact due to intrinsic 'wobbly 
swimming' rather than true (thermal) Brownian rotation |7D. 

An E. coli cell performs chemotaxis {i.e. swims up a chemoattractant gradi- 
ent) by alternately running and tumbling and biasing the length of the runs. The 
run durations are exponentially distributed with a rate constant ('tumble rate') of 
A w 1 s^', and a corresponding mean run duration of w 1 s [§]. Tumble dura- 
tions are also exponentially distributed but with a much shorter mean duration of 
w 0.1 s. An coli cell is too small to detect spatial differences in the concentra- 
tion of a chemoattractant on the scale of the cell length, so it performs temporal 
comparisons instead. It continually 'measures' the concentration of chemoattrac- 
tants (typically nutrients) in its environment such as serine and aspartate. If the 
concentration over the past second is higher than the concentration over the pre- 
vious three seconds, the tumble rate is reduced and the expected run length is 
extended. Thus, the cell performs a biased random walk and gradually drifts to- 
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ward regions of high chemoattractant concentration. The cell's mean velocity up 
the chemoattractant gradient is termed the 'drift velocity' . 

There have been a number of theoretical studies on run and tumble chemotaxis 
with no background flow |29, 14, 15, 16, 12, 25], and on the role of bacterial shape 
in chemotaxis 133, and references therein]. Some theoretical work has also been 
published on run and tumble chemotaxis with background flows. Bowen et al l ldl 
simulated bacterial chemotaxis toward a neutrally buoyant phytoplankton cell ex- 
uding dissolved organic carbon in an unsteady shear flow. Luchsinger et al ll25l 
simulated a similar situation, but also investigated the case in which bacteria re- 
verse direction rather than tumble. Bearon and Pedley [5] and Bearon 0] devel- 
oped an advection-diffusion equation describing run and tumble chemotaxis in a 
shear flow. The work of Bearon and Pedley |5] provided a particular motivation 
for the present study, since it showed that time-delays in the bacterial response can 
have a strong effect on drift velocity in a shear flow. 



The present work is in one sense an extension of Locsei l25[l to include shear 
flows, and in another sense an extension of Bearon and Pedley pj to include tem- 

poral comparisons, non-spherical cell shape, persistence, and Brownian rotation. 
Section|2]introduces the model and its basic assumptions, and section [3] describes 
in more detail the response function that is used to model the temporal compar- 
isons performed by the cell. Section |4] presents a general analytic framework for 
calculating the drift velocity. The remainder of the paper is devoted to specific 
calculations of the drift velocity while taking account of (i) temporal comparisons 
performed by the cell, (ii) persistence of direction, (iii) Brownian rotation, and 
(iv) non-spherical cell shape. The drift velocity under the combined effects (i), 
(ii) and (iii) is computed analytically in section |5l and the drift velocity under 
the combined effects (i) and (iv) is computed semi-analytically in section |6] The 
case where all effects (i)-(iv) are present is treated numerically by Monte Carlo 
simulation in section|7] Section[8]summarises the key findings. 

2 Outline of model 

Consider a cell performing run and tumble chemotaxis, swimming in an unbounded 
fluid with a chemoattractant concentration gradient Vc 1 1 z and with a background 
shear flow u = Jzx where x and z are unit vectors in the x and z directions and y 
is the shear rate. For convenience, we shall frequently parameterise the strength 
of the shear flow by 12 = y/2 rather than y. We assume that during a run the cell 
swims at constant speed Vs relative to the fluid. The cell's swimming direction 
is described by a unit vector e that changes with time due to tumbles, Brownian 
rotation, and the shear flow. 

During a run, the probability that the cell tumbles in the next time interval dt 
is k{t)dt, where A is the 'tumble rate'. We assume that the cell alters its tumble 
rate linearly in response the chemoattractant, according to 



where Ao = 1 s is the baseline tumble rate, and the fractional change in tumble 
rate is given by 




A(f)=Ao [1-4(0], 



(1) 




(2) 
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where c(f') is the chemoattractant concentration experienced by the cell at time 
t' , and R is the cell's 'response function'. R may be thought of as the impulse 
response of the tumble rate, since it describes the way that X changes when the 
cell is subject to a delta function impulse of chemoattractant concentration. The 
form of tumble modulation given by (O has been used in several earlier models 
iH, Q [H, Ell and it is motivated by the experimental results of [9] and |30]. 
The form of R is discussed in section [3j Our analysis will be restricted to 'weak 
chemotaxis', meaning small fractional changes in the tumble rate, i.e. \A{t) \ <^ 1. 
Physically, weak chemotaxis corresponds to a shallow chemoattractant gradient. 

Brownian rotation causes random re-orientation of the cell swimming direc- 
tion, so that in between tumbles the probability density function / of the swim- 
ming direction e evolves according to the Fokker-Planck equation 

|J+V,-[fi)(e)x*?/]=Z)RV2/, (3) 

where fi)(e) is the deterministic angular velocity of the cell due to the shear 
flow, Z)r is the rotational diffusion coefficient and Vj, is the gradient operator 
in direction space. The form of (0{e) depends on what assumptions are made 
about the cell shape. As discussed in the introduction, the rotational diffusivity 
of the cell may be due to intrinsic randomness in its swimming motion, rather 
than true Brownian rotation. Throughout this paper we shall set either Dr = or 
Z)r = 0.062 radians^s"^, where the latter value is consistent with what has been 
measured in experiments on E. coli [6]. 

When a cell tumbles, its choice of new direction is governed by a probability 
distribution which is axisymmetric about the initial direction. We allow for the 
tumbles to exhibit directional persistence, so that the expected scalar product of 
the swimming directions e{0~) and e{0^) immediately before and after a tumble 
at time f = is given by 

E[e(0-)-e{0+)] = ap, (4) 

where E denotes an expectation value and Op is the 'persistence parameter' . Ex- 
perimentally, 0(p w 0.33 1 6]. We treat tumbles as instantaneous and neglect any 
rotation of the cell caused by the shear flow during the tumble. 

Our aim is to calculate the chemotactic drift velocity. Let Zf be the z location 
of a cell at the end of a run, relative to its position at the beginning of the run, 
and let tf be the duration of a run. Since we treat tumbles as instantaneous, the 
chemotactic drift velocity Vd is 

Vd = £[zf]/£[ff]. (5) 

While our model assumes an unbounded domain, we note that for a cell swimming 
in a bounded domain of length ^ Vs/Ao, provides a measure of the transient 
velocity up the chemoattractant gradient before the cell encounters the boundaries. 



3 Response function 

In writing (O, we assume that the tumble rate is a linear functional of the con- 
centration history seen by the cell. The validity of this assumption has never been 
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Fig. 1 The response function used in the results section of this paper, originally proposed by 
Clark and Grant UJ]. 



directly assessed, and the chemotactic response function R has never been directly 
measured. However, experiments by Block et al and Segall et al |31] indicate 
that the bias of a single flagellar motor (the probability of counter-clockwise ro- 
tation) is a linear functional of the concentration history, thus providing indirect 
evidence for ([T|. These authors monitored the direction of rotation of a single ro- 
tary motor on an E. coli cell while delivering small impulses of chemoattractant to 
the fluid around the cell. Repeating the experiment multiple times and with differ- 
ent cells, the experimenters determined the motor bias as a function of time (see 
for example figure 1 of Segall et al 13 1.1 ). The impulse response of the motor bias 
is double-lobed, with the bias rising above the baseline for the first wis after 
the delivery of the impulse, falling below the baseline for the following w 3 s, 
and then returning to baseline. The two lobes of the response have equal area. 
Furthermore, it was found that the responses to other time-series of stimuli {e.g. 
ramp or sinusoidal changes in chemoattractant concentration) are consistent with 
the cell behaving as a linear system, so that the motor bias is well described by 
the convolution integral of the stimulus with the impulse response {c.f. equation 
[T|). The primary exception to this linear behaviour is that for small changes in 
chemoattractant concentrations, cells respond to increases in concentration but 
not to decreases; we neglect this nonlinearity in our analysis. 

Following the lead of previous papers I29l[l4l [l2il . we adopt the view that in 
the absence of experimental measurements of R, a convenient assumption is that R 
has a similar shape to the impulse response of the individual motor bias reported 
in Block et al |9] and Segall et al [31]. This means that (i) R should be composed 
of a positive lobe followed by a negative lobe of equal area so that 

/ R{t)dt = 0, (6) 

JO 

and (ii) R{t) should decay to zero for t greater than about 4 s. An additional con- 
straint is that R must be small enough for our linear analysis to be valid. Certainly, 
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R must be small enough that A never becomes negative, which necessitates 



|4ff)L.x=vJVc| 



/ tR{t)dt 
Jo 



< 1- (7) 



The analysis that will be presented in this paper is asymptotically accurate as 
We shall use the response function 

m = 3^e-^'[l - Aot/2- (Aot/2)2], (8) 

with e ^ 1 in order for our linear analysis to be valid. This R, plotted in figure[Tl is 
a theoretically motivated response function proposed by Clark and Grant lfl2ll that 
has the above-mentioned properties (i) and (ii), and it matches the experimentally 
measured motor bias reasonably well. It may seem odd that |Vc| appears on the 
right hand side of ([Hi, since 7? is a property of the cell rather than its environment. 
However the factor of l/|Vc| in R simply reflects the fact that we have used e 
to parameterise the strength of the combined effect of the cell's response R and 
the chemoattractant gradient Vc on 4 (/ ) . Note that in the computations we do not 
have to specify a value for | Vc|, since R and | Vc| appear together as a product in 
the equation for A{t). 



4 Partial calculation of drift velocity 

In this section we present the first part of the analytic drift velocity calculation, 
which is common to both the spherical and non-spherical cell cases. Calculating 
Vd is non-trivial because of the interdependence of the tumble rate and the path 
taken by the cell. The tumble rate at any time depends in principle on the entire 
path history of the cell through (O, while the path of the cell depends in turn on 
the tumble rate. The method of solution is the same as that used in Locsei |25], 
so we shall omit some details. As noted earlier, in order to make the analysis 
tractable, our analysis will be restricted to weak chemotaxis, i.e. <\A{t)\<^ I. 
In this case, the expectation values of the run displacement Z{ and duration tf satisfy 
E[zf] = VsAoO(4) and E[tf] = [1 + C»(zi)]/Ao, so, from ©, 

Vd = Ao£[zf]+v,0(zi2). (9) 

We shall neglect terms that are 0{A^). 

Consider a run commencing at time / = at location z = 0. During the run, the 
cell swims in a random walk governed by rotational Brownian motion until the 
run terminates with a tumble at time tf. Note that the expected stopping location 
of a terminated random walk with a stopping rate Astop is the same as the expected 
first event location on an unterminated walk with an event rate Aevent = ^stop ■ Thus, 
in calculating i?[zf], it is permissible to treat tumbles for f > as events that have 
no effect on the cell's motion, and treat the tumble at f = ff, z = Zf as a first event 
(the first tumble in ? > 0). The utility of this treatment is that we may conceptually 
break the expectation E in E[zi\ into two consecutive operations. First, assuming 
a given path z{t) : —°° <t<oo taken by the cell, one calculates the conditional 
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expectation of Zf for that path. Second, one takes the expectation over all such 
paths to obtain E [zf] , with the understanding that in the t <0 section of a path the 
cell is subject to reorientations due to both Brownian motion and tumbles, whereas 
in the t > section of a path the cell is subject to reorientation due to Brownian 
motion alone. Writing out the two expectations in symbolic notation. 



E[zf] — -Epaths 



dtZpath{t)Ppath{t) 



(10) 



where the £paths denotes an expectation over paths, Zpath(f) denotes the position 
of the cell at time f on a particular path, and /?path(0 is the probability density 
function for the tumble time ff on a particular path. Since tumbles for f > are 
treated as having no effect on cell motion, paths are independent of t{ and one is 
free to take the path expectation inside the integral over tumble times. Substituting 
(llOb into (|9]l and dropping the 'path' subscript for brevity, one then has 



(11) 



Vd = Ao / dtE[z{t)p{t)]. 
Jo 

The probability density function p(t) for the tumble time is given by 

p(f)=A(f)exp -/ k{t')dt' 
Jo 



(12) 



where A(f) is the path-dependent tumble rate at time t, given by (dJ. 

We shall calculate the drift velocity for the case where R is given by a Dirac 
delta function, 

R{t)=Ad{t~T), (13) 

and later generalise to an arbitrary response function. Expanding (II lb in powers 
of A, keeping only the linear term, then using equation ( [13]) for R and rearranging, 
one finds 

Vd = Ao/' dts-^'E[w{t)] + ?i^A\Wc\ [ dt' [ dts~^'E[w{t)z{t' -T)], (U) 
Jo Jo Jt' 

where w(f) =t/z(f)/t/f. Writing as an integral of w(f), writing w(f) as Vse(t).z, 
( fT4] i becomes 



Vd ■ 



Ao / dte-^'z-E[eit)] 
Jo 

+?iiA\Vc\ [ dt' ( dt [ dt"e-^'z-E[e{t)e{t")]-z, (15) 
JO Jt' Jo 



Equation (115b contains the unknown terms E[e{t)] and E[e{t)e{t")]. The cal- 
culation of these terms under different sets of assumptions is the subject of the 
following two sections. 
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5 Analytic calculation of for a spherical organism 

5.1 Derivation 

For the case of a spherical organism, one can solve for Vd analytically, taking full 
account of temporal comparisons, persistence, and Brownian rotation. Once again, 
the method of solution is similar to that in Locsei lisll . so just an outline of the 
derivation is sketched here. 

First consider the case with no shear. In this case, the expected swimming 
direction for times f > is 

i?[e(f)]=e-2^R'i?[e(0+)], (16) 



where e(0+) is the cell's swimming direction at the beginning of the run 12511 . The 
decaying exponential reflects the random re-orientation of the cell by Brownian 
rotation. The time autocorrelation function for swimming direction is 



g2Z)R(ri-r2)2 ifO<fi <f2 



E[e{t2)e{H)] = | BieP^R+^H-^D'-e-^^R'^l if fi < < f2 ^^^^ 

where 1 is the identity matrix. Note that the second case in dlTl l has a factor Op 
due to the tumble at f = 0, an exponential decay term featuring Dr and Xq due 
to tumbles and Brownian rotation in the interval fi < f < 0, and an exponential 
decay term featuring just Dr due to Brownian rotation in the interval < f < f2- 
For Reynolds numbers Re <C 1, a spherical object in a shear flow rotates with an 
angular velocity Q) equal to half the vorticity i23ll . When shear is included, one 
can show that ( 1161) and ( 117b generalise to 

£[e(0]=e-2DRt^?(t).E[e(0+)], (18) 

and 

zrr / N / M f ie2°R('l-'2)7?(f2-fi) if0<fl<f2 

E[e{,t2)e[,h)\ I 5ie{2^R+^[i-«pl)'.e-20R'2^(f2 if < < f2 ^ ^ 

where R{t) is the rotation matrix corresponding to a rotation by angle \(0\t about 
the unit vector 6). 

It remains to find an expression for £[^(0^)]. Since the swimming direction 
after a tumble is chosen from an axisymmetric distribution around the swimming 
direction before the tumble, @ gives 

E[e{0+)] = apE[e{0-)]. (20) 

Furthermore, if the cell has been swimming for a time fswim ^ 1 /■^, then the ex- 
pected swimming velocity just before a tumble does not change from one tumble 
to the next, so 

E[e{0+)] = a^E[e{tf)] 

= Op f dtE[e{t)p{t)]. (21) 
Jo 
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After expanding p{t) in powers of A and performing some simplification, one 
finds ^ 

E[eiO+)]=vM^c\A(B-—l) -{h-h) (22) 



where 



Op 

/■oo f[t-T) 

dt dt'Q-^'E[e{t)e{t% (23) 

./o Jo 

h = ^l dt' [ dt f ^ dt"e'^'E[e{t)e{t")], (24) 
Jo Jt' Jo 



and 



B = ?io[ t/fe-(^+2^R>/?(f). (25) 
./o 

Now we have expressions for E[e{t)] and E[e{t)e{t")], and performing the 
necessary back-substitutions into ( 1151) yields 

Vd =^^sphere(r), (26) 

where 

, ^ .-n^^R+^)|Vc|v^(«p-l)Aoj.i 

/CsphereUJ 3 [(2Z)r + Ao)2 + [(2Z)r - OpAfl + Ao)2 + i22]2 ' ^^'^ 

with 

^gr(2DR+Ao)^;^^^2_grapAo[(2Z)R + Ao)2 + I22]/,3 (28) 

= (2Z)r + Ao) (2Dr - OpAo + Ao)2 + [(2ap - 3)Ao - 6Z)R]i22 (29) 

b3= [(2Z)R-apAo + Ao)^-I2^]cosI2r 

-2(2DR-apAo + Ao)i2sinX2r. (30) 

The drift velocity for a general response function is then 

Vd = j^ dTR{T)k,p^,,,{T). (31) 



5.2 Consistency with known results 

Equations jZTl) and ([3T| are consistent with existing literature in the appropriate 
limits. Locsei ll25h considered chemotaxis in the absence of shear. When £2=0, 
(|27] | simplifies to 

v2|Vc|V-(^+2DR)r(i_c^)[(^^2DR)eW-Aoape(^+2DR)r] 

^ ^ 3(Ao + 2DR)[^(l-ap) + 2DR]2 

(32) 

consistent with f2y]. Bearon and Pedley pj calculated the chemotactic drift ve- 
locity for a spherical cell swimming in a uniform shear, under the assumption that 
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the cell detects the chemoattractant gradient directly (rather than through temporal 
comparisons) and modifies its tumble rate according to 

A =Ao(l-e<?-z), (33) 

with |e| <^ 1. We note that within our framework ( [33] l is equivalent to setting 

R(t)- ' 5{t-T)-5{t-T-At) 

^^'^-^ At ' ^^^^ 

and taking the simultaneous limits At ^ and T ^ 0. Substituting 041) into d^TT) 
and taking the aforementioned limits, one finds the drift velocity for a gradient- 
detecting spherical organism: 

^ £VsAo(ap-l)[Ao(ap-l)-2£)R] 

3[(Ao(ap-l)-2DR)2 + i22]. ^^"^ 

Setting Op = and = 0, one finds 

''=3{?.§+n^y ^^^^ 

consistent with equation 30 of Bearon and Pedley 



5.3 Results: Vd for a spherical organism 

Using response function ([8]l, the drift velocity for a spherical organism in 
flow is 

Vd _ {ap~l)?ii{b5[2DR-{ap-2)?^]n^-b4 + be} 

where 

b4 = [4£>R + (5 - 2ap)Ao] [2Z)r - (Op - 2)Ao]3 (2Z)r - OpAo + Ao), (38) 

bs = 16D| + 4(ll-4ap)AoZ)R + (4ap2-22ap + 25)Ao^ (39) 

^6 = 3[4Z)R + (3-2ap)Ao]X24, (40) 

and 

b^ = 9{[(ap - 2)Ao - 2Dr]2 + Q^}^[{2Dk - OpAo + Ao)^ + Q^] . (41) 

Figure |2] shows the drift velocity Vd as a function of the shear strength Q for a 
variety of persistence parameters Up. As one might expect, the shear flow reduces 
the drift velocity. More surprisingly, the drift velocity is negative for sufficiently 
strong shear, so that the cell swims down the chemoattractant gradient rather than 
up it. Also, the larger the persistence parameter, the more susceptible the cell is 
to shear. For instance, with a persistence parameter of Op = — 1, Vd becomes neg- 
ative for ^2 beyond w 1.5 s~^, whereas for Op = 0.78, the drift velocity becomes 
negative for Q beyond w 0.4 s~^. Locsei [25J calculated that in the absence of 
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Fig. 2 Drift velocity versus shear strength, for persistence parameter Kp = — 1 (solid line), Op = 
0.33 (dashed), and Bp = 0.78 (dotted), with Ao = 1 s"' and Dr = 0.062 radians^s"'. Note ttp = 
0.33 is the measured value for E. coli j^. 



shear, Op w 0.78 is the value of persistence that maximises the drift velocity. In 
that paper, it was speculated that the measured persistence of Ofp w 0.33 < 0.78 
in real E. coli \6\ might reflect a compromise between transient and steady state 
performance. The new results here suggest another reason why < 0.78 for real 
E. coli: there is a compromise between performance in a stationary fluid and per- 
formance in a sheared fluid. Intuitively, the disadvantage of a high persistence 
parameter when performing chemotaxis in a sheared fluid arises because when 
a cell is rotated by the shear to face down the chemoattractant gradient, a single 
tumble is insufficient to reorient it to face up the gradient. 



6 Semi-analytic calculation of Vd for an elongated organism 

Whereas a spherical particle in a shear flow rotates with constant angular veloc- 
ity, a non-spherical particle rotates with an angular velocity that depends on its 
instantaneous orientation. This complicates the calculation of drift velocity, but 
nonetheless one can still find an analytic integral expression for the drift velocity 
of a non-spherical cell in the absence of Brownian rotation and persistence, using 
the additional simplification of treating the cell as a prolate spheroid. In principle, 
one could include Brownian rotation and persistence in the calculation, but this 
would involve solving the full, time-dependent Fokker Plank equation, and we 
shall not attempt this here. 



6.1 Modelling an E. coli cell as a prolate spheroid 

In this section we show that an E. coli cell in a shear flow rotates in approximately 
the same manner as a prolate spheroid (ellipsoidal rod) of slenderness ratio T] w 9, 
where 77 is defined as the ratio of major to minor axis. 
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We begin by considering the motion of a rigid body in a linear velocity field. 
A general, undisturbed linear velocity field a (x) has the form 

u{x) =U + Q xx+E-x, (42) 

where [/ is a uniform velocity, Q is an angular velocity equal to half the vorticity, 
and £ is a second rank symmetric strain rate tensor. The shear flow considered in 
this paper is a special case of a linear flow, with 

U = 0, 

Q = (0,12,0), 

f on\ 

£= 0, (43) 
\Q0 J 

where Q = y/2. 

A swimming E. coli cell has a Reynolds number of Re w 10~^ so the 
hydrodynamics of its motion are well described by Stokes flow. By the linearity 
properties of Stokes flow, the instantaneous angular velocity 03 of any rigid body 
in a linear flow takes the form 

(Oi = ni+BijkEjk (44) 

(using the Einstein summation convention), where 5,^^. is a dimensionless third 
rank constant tensor which is symmetric in its second and third indices, and the i, 
j, and k components refer to axes fixed in the body. For a sphere, B = 0, reflec ting 
the fact that a sphere does not rotate in a pure straining flow. Bretherton illh 
showed that for a general axisymmetric body 5,,^ has only a single free parameter. 
Taking the z axis to be aligned with the body's symmetry axis, an axisymmetric 
body has 

-^123 = -5i32 = S2I3 = ^231 = a, 

Bijk = for all other ij,k. (45) 

Bretherton i fTTI) further showed that, provided |o;| < 1 (which is true for most 
physically realistic bodies), a body with Bjjk of the form ( 1451) rotates in a shear 
flow in precisely the same way as a spheroid of slenderness ratio 



77 = V(l + a)/(l-a)- (46) 

This result is useful because the motion of a prolate spheroid in a shear flow was 
analysed by Jeffery (tS], who showed that the motion consists of a closed orbit of 
simple analytical form, now commonly referred to as a 'Jeffery orbit' . 

We shall consider a model E. coli cell composed of a spherical body rigidly 
attached to a helical flagellar bundle, modelled as a solid body of uniform circular 
cross-section, and show that despite the model cell being non-axisymmetric its 
Bjjji tensor is almost of the form ( |45]) . Our model cell consists of a sphere of radius 
a centered at the origin with respect to axes fixed in the cell and a left-handed helix 
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of radius r, pitch (helical wavelength) h and length L, such that the position jc at a 
point on the helix is given parametrically by 

x{x) = r{nico^X - ni^mx) +n3 (|« + ^^ : (47) 

where ;^ e [0,2nL/h] and«i,«2,«3 are a right-handed set of unit vectors fixed in 
the cell such that n^, points along the helix axis away from the sphere. 

It is sufficient to treat the helix and sphere as being fixed in position relative 
to one another, so that the cell is rigid rather than swimming. By the linearity 
property of Stokes flow, the angular velocity of a swimming cell in a shear flow 
is equal to the angular velocity of a swimming cell in a stationary fluid, plus the 
angular velocity of a rigid cell in a shear flow. Since D.i= for a swimming 

cell in a stationary fluid (because the cell swims in an approximately straight line), 
Q.\ and are the same for a swimming cell in a shear flow as for a rigid cell in a 
shear flow. 123 is not relevant to calculations of drift velocity, since rotation about 
«3 does not alter the swimming direction. 

For a given background flow, the translational velocity V and angular velocity 
0) of the body are found by setting the net force and net torque on the body to zero. 
We neglect the interaction of flows around the helix and the sphere. For a given 
motion of the cell through a given background velocity field, we assume that the 
total viscous force on the cell is simply the force on the sphere in the absence of 
the helix plus the force on the helix in the absence of the sphere, and we simi- 
larly assume simple addition of the viscous torques. This is not an unreasonable 
approximation since the body diameter (l-2jUm) is a small fraction of flagellar 
bundle length (5-lOjUm). To calculate the force and torque on the sphere, we use 
the well known exact expressions [see for instance i2j1. To calculate the force and 
torque on the helix we use resistive force theory flSll . which assumes that the 
components of the viscous force on an element of the helix are proportional to the 
same components of the fluid's velocity relative to that element, but with differ- 
ent coefficients of proportionality for normal and tangential components. Various 
pairs of coefficients have been proposed for a simple flagellum of uniform circu- 
lar cross-section, and we use the more accurate coefficients due to Lighthill |24l 
rather than the original coefficients proposed by Gray and Hancock jlSll . These 
coefficients have as parameters the helix pitch h and the cross-sectional radius b 
of the helix bundle. Resistive force theory is a crude approximation of true flag- 
ellar hydrodynamics, since it ignores interactions between neighbouring elements 
of the flagellum, but it has nonetheless been shown to give reasonable results when 
applied to swimming spermatozoa, for example |24]. 

Given the above method for calculating the motion of the body in a given flow, 
one can determine the Biji^ tensor from ( 1441) by setting £2=0 and considering 
the angular velocity CD that results from each of the following elementary strain 
matrices, in which all Ejj except those mentioned are zero: En = 1; E22 = 1; 
£33 = 1; £12 = E21 = 1; E23 = -E32 = 1; £^31 = £^13 = 1- 

In principle one can determine S/yj. in symbolic form with a, b, r, L, and h as 
parameters, but in practice the algebra required to perform this task is unwieldy 
and so one must choose numerical values for the parameters before calculating 
Biji^. A typical E. coli cell body is rod shaped with a diameter of w 1/im and a 
length of w 2/im, and possesses about 6 flagella [6]. The length of each flagellum 
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is typically in the range 5-10/im, and the diameter of each flagellar filament is 
w 20nm [28, pp. 70-83]. The pitch of the helix is w 2.5jUm and the helix radius is 
w O.lSjim |32]. On the basis of these measurements, we take the dimensions of 
our model cell to be a = OJjim, r = 0.25;Um, L = Sjim, h = 2.5/im, b = \f(j x 
lOnm. With these dimensions, calculation of the reveals that, for the model 
cell, they come close to possessing the properties expected for an axisymmetric 
body [eq. (l45Tl1. in that 

~ 5i23 = — 5i32 ~ ^213 = S231 ~ 0.974, 
S333 » 0.501, 
I < 0.05 for all other «, j, k. (48) 

The most notable breaking of symmetry is that S333 is not small (whereas it is zero 
for an axisymmetric body) but this only affects Q.^, and thus has no bearing on the 
swimming direction. All other elements of that are zero for an axisymmetric 
body have magnitude < 0.05 for the model cell (and all but one, S133 = —0.0445, 
have magnitude < 0.01), so discrepancies between the angular velocities of the 
model cell and an 'equivalent spheroid' are of order 5%. Using ( 146b , the slen- 
derness ratio of the 'equivalent spheroid' that rotates in a similar manner to the 
model cell is 77 w 9. For other plausible values of a and L, 77 ranges from T] w 4 
(for a = ljum, L = 5jUm) to 77 w 15 (for a = 0.5jUm, L = 10/Jm). 



6.2 Integral expression for drift velocity of prolate spheroid 

Neglecting Brownian rotation and persistence, one can find an analytic integral 
expression for the drift velocity. The direction faced at the beginning of each run is 
taken from the uniform distribution on the unit sphere, so by symmetry ii[e(f)] = 
and the first term in ( 1151 ) vanishes. Also, since any correlation in swimming 
direction is destroyed by the tumble, E[e{t)e{t")] = for f" < < t, which allows 
some simplification of the second term in ( 1151) . One finds 

Q24|y I ^00 r(t-T) 

^ ' / dt^-^' dt'{t-t'-T)E[w{t)w{t')]. (49) 

vj .It Jo 

It remains to find an expression for the velocity autocorrelation function E[w{t)w{t') 

Since it was established in section|6T|that an E. coli cell rotates in a shear flow 
with very nearly the same angular velocity as a prolate spheroid, we henceforth 
treat the cell as a prolate spheroid that swims along its long axis. Define direction 
angles and (j) such that 

e-z = &m9 cos (/) (50) 
e Jc = sin0 sin0 (51) 
ey = cosd. (52) 

In a shear flow, neglecting Brownian rotation, the cell follows Jeffery orbits, and 
the angles 6 and (/) evolve according to ll22h 

tan (/) = 77 tan(af - j!3o) , (53) 
,2 



,2^ r 

JCq ( 7] 2 cos2 (/) + sin"^ (/) ) 



tan^ e = ; , , . 2 , . : (54) 
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where Ko and j3o are constants of integration depending on the initial orientation, 
and 

2X2 T? 

In terms of the initial orientation angles do and at time t = 0, the constants of 
integration are 

ft ( </'o ) = - arctan( T] cos 00 , sin ) , (56) 

2 

tan"^ 0o(7]"^cos"^(/)o + sm 00 ) 

where the arctan function is treated as taking two arguments so as to return an 
answer in the correct quadrant. For instance, if i// = arctan(ji:i ,ji:2), then cosy/ = 
xi{x1+x]y^^^ and sin^ = X2{x\+xl)~^^^. 

Let w{t;eo) be the cell's velocity projected along the z axis at a time t > 0, 
given that the cell's swimming direction at the beginning of the run was cq. Ma- 
nipulation of (|50ll, (l53]l and ( l54l) gives 



w{t;eo) = Vssin0cos0 

VsCos(a/ — j3o) 



[ka +cos2(af - Po) + 772 sin2((7f - ft)] 



1/2- 



(58) 



E[w{ti)w{t2)] may be written as an explicit average over all possible initial run 
directions; 

J i-n 1-271 

E[w(ti)w{t2)] = — c/QosinSo / ii(/>ow(fi;<?o)w(f2;<?o), (59) 
471 Jo Jo 



where the integrand is treated as a function of 6q and ^q. Using (156b . we can 
transform this last equation into to 

47rJo Jo cos2j3o + ?7"^sm ft 

where the integrand is now treated as a function of do and ft, and there is an 
accompanying expression for Kq in terms of do and ft instead of Oq and 0o- 

^fo R\ cos2/3o + T]2sin^ft 

K{)(0o,ft = • (61) 

tan'^ t^o 

Equations ( |58] l. ( |60] l. and (|6T| now give a complete expression for E[w{ti)w {12)]. 
The right hand side of ( 1601 ) cannot be integrated analytically for general 77, but can 
be evaluated by numerical quadrature. In the limit of weak shear, one can expand 
the integrand in powers of X2 /Ao, taking f 1 and t2 to be of order 1 /Aq, and integrate 
to find 

E[w{H)w{t2)] = f-j^^0^{-2i3r^' + Sil' + 24)Ht2 

+ (12774 + 1 1772 + 12)(ff + + 0{Q^/?^). (62) 
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Note that ( |62] | contains only even powers of Q. since by symmetry the drift velocity 
does not depend on the sign of Q. . 

For a general response function, the drift velocity is given by 



Vd 

where 



dTR{T)k,^k.roi<i{T), (63) 



Q2|y„| ^oo r{t-T) 

fcsphero.d(2^) = dtc-^' dt' {t-t'- T)E[w{t)w{t% (64) 

withii[w(f)w(f')] given by ( |60] |. 



6.3 Consistency with known results 

de Gennes calculated the chemotactic drift velocity in the absence of a shear 
flow, and without Brownian rotation or persistence. Setting X2 = 0, the velocity 
autocorrelation function is simply i?[w(fi),w(f2)] = Vs/3, and ( |64] | and (|63]) reduce 
to 

Vd = ^ rR{t)c-'o'dt, (65) 

consistent with equation 16 of de Gennes fl?]. As mentioned in section 15.21 
Bearon and Pedley |5] calculated the chemotactic drift velocity for a spherical 
cell swimming in a uniform shear, under the assumption that the cell detects the 
chemoattractant gradient and modifies its tumble rate according to ( 1331 ). We can 
treat this case as follows. For a sphere, T] = 1, and the velocity autocorrelation 
function simplifies to 

E[w{ti)Mt2)] = jCos[n{ti-t2)], (66) 



and substitution of S66i into dMI) yields 

c- [{Q^ -^^) cos QT + 2Q sin QT] 
'"^^'^ 3(A2 + I22)2 ■ ^^^^ 

As in section 15 .2[ gradient detection is modelled by using the response function 
( 134b and taking the appropriate limits, and one finds that once again the drift ve- 
locity reduces to ( l36b . consistent with 



6.4 Results: Vd for an elongated organism 

For weak shear {Q <^ Ao), we can substitute ( |62] | and ^ into ( |64] | and ( 1631) to find 

Vd 



= e 
evs 



5 , , fn^ 



(68) 
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Fig. 3 Drift velocity as a function of the shear strength parameter 12 for equivalent slender- 
ness ratios rj = 1 (solid line), r) = 4 (dashed line), t] = 9 (dotted line), and r) = 15 (dot-dashed 
line). 



where 

, , 243T]4 + 15l772-44 

Note that ^(77) is bounded and monotonia ally increasing for ?7 > 0, with g{°°) — 
81/140. In fact, g{4) w 0.92g(oo), indicating that for 12 < 1 and T] > 4, Vd is 
approximately independent of T] . 

For general Q and 77, no analytic simplifications are possible and we evaluate 
the integral ( 163b numerically using the 'Cuba' integration library developed by 
Hahn [20]. Figure |3] shows the drift velocity as a function of the shear strength 
Q for a range of values of the slenderness ratio rj . Note that for 12 < 0.5 and r/ > 4, 
Vd is approximately independent of T], consistent with the analysis in the previous 
paragraph. Also, for any given T], and for sufficiently large Q, the cell exhibits 
'negative chemotaxis' i.e. Vd is negative. For instance, if T] = 9 then Vd < for 
Q > 3. The minimum value of Q for which Vd < depends on T]; for larger T], Vd 
remains positive up to a larger value of Q . 

This last finding can be explained as follows. A key feature of the Jeffery orbit 
of a prolate spheroid with large slenderness ratio is that that there are lengthy 
'pauses' with the spheroid's symmetry axis almost parallel to the streamlines, but 
the orientation reverses periodically. From ( |53l l and ( |55l l we see that the period 
Tjef of a Jeffery orbit is 

7}ef. = 7^(772 + l)/(r2T]), (70) 

which is a monotonically increasing function of t] for 77 > 1. Let i2crit be the low- 
est value of Q for which Vd = 0. In order for Vd to be less than or equal to zero, the 
sign of w must change during a typical run, and thus Tjef. must be comparable to 
or smaller than the mean run duration I/Aq. Thus, we expect that the dependence 
of i2crit on T] scales roughly as 

i2crit/Ao-?(T]) = (T]2 + l)/T7. (71) 

In fact, we find that (i2crit/Ao)/^(77) = 0.42,0.35,0.33,0.32 for 77 = 1,4,8,16 
respectively, indicating that the scaling ( 171b is approximately correct. 
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Fig. 4 Examples of agreement between analytic theory and Monte Carlo simulation re- 
sults. Solid line shows analytic result for a spherical organism with Op = 0.33, Dr = 
0.062 radians^s^', e = 0.3, and Ao = 1 s^'; solid squares are Monte Carlo results for the same 
set of parameters. Dashed line shows analytic result for an elongated organism with 7] = 9, 
(Xp =0, Dr = radians^s^', e = 0.3, and Ao = 1 s^'; empty squares show Monte Carlo results 
for the same set of parameters. Error bars show the standard error of the mean for Monte Carlo 
results. 



7 Calculation of by Monte Carlo simulation 

The purpose of using Monte Carlo simulation here is twofold. First, it gives a 
means of checking the analytic results derived in earlier sections. Second, it al- 
lows us to calculate in the regime which was not analytically tractable, namely 
the regime including all four effects of (i) temporal comparisons, (ii) persistence, 
(iii) Brownian rotation, and (iv) elongated cell shape. One can then check that 
the conclusions drawn from the analytic calculations are still valid when all four 
effects are taken into account. 



7.1 Monte Carlo algorithm and test cases 

The algorithm for calculating Vd by Monte Carlo simulation is straightforward. For 
each set of parameter values, 100 simulations are run, each of duration 10^/Ao- 
Time is discretised into steps of size At = 0.01 /Ao. At each time step, the tumble 
probability ([T| is evaluated by numerical quadrature, and a random number gen- 
erator is used to decide whether the cell tumbles. If the cell tumbles, then the new 
direction is chosen from an axisymmetric distribution about the old direction, such 
that the new direction makes an angle of arccos(ap) with the old direction. If the 
cell doesn't tumble, then Brownian rotation is simulated by giving the cell a new 
direction chosen from an axisymmetric distribution about the old direction, such 
that the new direction makes an angle of arccos(l — ID^^At) w I^/DrAI with the 
old direction, and in addition the cell's position and direction of motion is incre- 
mented according to the Jeffery orbit equations (23}. A computer code was written 
to implement the above algorithm, and its output was compared to the analytic re- 
sults derived in earlier sections. As show in figured the analytic and simulation 
results agree well. 
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Fig. 5 Monte Carlo results for drift velocity Vd as a function of shear strength Q, for persistence 
parameters of 0(p = — 1 (solid line), Op = 0.33 (dashed), and Op = 0.78 (dotted). In each case, 
T] = 9, Z)r = 0.062 radians^s^', Ao = 1 s^' and £ = 0.3. Standard errors (not shown) are of the 
same magnitude as those in figure|4] 



7.2 The effect of persistence on drift velocity 

In section 15.31 it was shown that for a spherical organism, a higher persistence 
OCp results in a greater drift velocity in the absence of shear, but a lower drift 
velocity in higher values of shear (figure |2]l. Using the Monte Carlo code, the 
same calculation was repeated, but for an elongated organism (t] = 9). As shown 
in figure |5] the results remain qualitatively similar. For instance. Op = 0.78 gives 
a higher Vd than Op = 0.33 for shear rates up to 12 w 0.2 s^', but for shear rates 
beyond that the situation is reversed. The main difference between figures |2] and 
|5] is that when elongated cell shape is taken into account, Vd never becomes as 
strongly negative as it does in the case of a spherical organism. 



7.3 The effect of organism shape on drift velocity 

In section |63] it was shown that a more elongated organism (larger rj) is able to 
perform chemotaxis in stronger shear flows than a less elongated organism (figure 
O, in the case where persistence and Brownian rotation are neglected. Using the 
Monte Carlo code, the same calculation was repeated, but including persistence 
(OCp — 0.33) and Brownian rotation (Dr = 0.062 radians^s"^). As shown in figure 
|6l the results remain qualitatively unchanged. For instance, for 77 = 1, Vd remains 
positive for shear rates up to 12 w 0.8 s"^ whereas for T] = 9, Vd remains positive 
for shear rates up to i2 w 1 .5s~ ^ . However, the effect of rj is less pronounced than 
in the case where persistence and Brownian rotation are neglected. The physical 
reason for this is as follows. In the absence of Brownian rotation, an elongated 
cell (t] ^ 1) has a lengthy 'pause' when it is almost parallel to the streamlines 
of the shear flow. The larger the value of T], the longer the pause, and the more 
robust Vd is in the face of shear. However Brownian rotation 'jostles' the cell out 
of alignment with the streamlines, puts a limit on the duration of the pause, and 
consequently reduces the robustness of Vd in the face of shear. 
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Fig. 6 Monte Carlo results for drift velocity as a function of shear strength Q, for slender- 
ness ratios of 77 = 1 (solid line), rj = 4 (dashed), and rj = 9 (dotted). In each case, Op = 0.33, 
Dr = 0.062 radians^s^', Ao = 1 s^' and e = 0.3. Standard errors (not shown) are of the same 
magnitude as those in figure|4] Results for 7] = 15 are statistically indistinguishable from results 
for T] = 9 (data not shown). 



8 Discussion and Conclusion 

In this paper we have used a combined analytic and numeric approach to calculate 
the drift velocity of a cell performing run and tumble chemotaxis in a uniform 
shear flow, taking account of (i) temporal comparisons performed by the cell, (ii) 
persistence of direction, (iii) Brownian rotation, and (iv) cell geometry. Our key 
findings are that (a) shear reduces the drift velocity, and a sufficiently strong shear 
may cause the cell to perform negative chemotaxis and swim down the chemoat- 
tractant gradient rather than up it, (b) in terms of maximising drift velocity, a high 
persistence parameter is advantageous in a quiescent fluid but disadvantageous in 
a shear flow, and (c) a more elongated body shape is advantageous in performing 
chemotaxis in a strong shear flow. 

Temporal comparisons of chemoattractant concentration performed by the cell 
constitute a form of time delay in its response. The cell biases its tumble rate ac- 
cording to the history of chemoattractant concentration it has experienced, rather 
than according to the instantaneous concentration or concentration gradient. Neg- 
ative chemotaxis arises from the combination of this time delay with rotation by 
the shear flow. For instance, a cell may commence swimming up the chemoat- 
tractant gradient on a particular run, but be rotated by the shear flow to swim 
down the gradient. Due to the time-delay in its response, the cell responds as if 
it were still swimming up the gradient, so it extends the length of the run, even 
though the net displacement is down the gradient. Conversely, a run with a net 
displacement up the chemoattractant gradient may be foreshortened. Lengthened 
runs down the gradient and foreshortened runs up the gradient lead to negative 
chemotaxis. Note that negative chemotaxis does not occur for an organism that 
responds to the instantaneous chemoattractant gradient; the drift velocity given by 
( 135b is non-negative. 

The possibility of negative chemotaxis was already derived by Bearon and 
Pedley |5], although the authors did not comment on it explicitly. They calculated 
the drift velocity for a spherical cell performing run and tumble chemotaxis in a 
shear flow, and investigated a form of time delay in the cell's response. In their 
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model, a cell could bias its tumble rate according to either the present angle be- 
tween the swimming direction and the chemoattractant gradient (no delay), or the 
original angle between the swimming direction and the chemoattractant gradient 
at the beginning of the run (delay). The authors found that in the delay case the 
drift velocity was significantly more reduced by shear than in the no-delay case, 
and equation 27 of their paper shows that in the delay case the drift velocity be- 
comes negative for sufficiently large shear. 

The present study shows that in terms of maximising drift velocity, a high 
persistence parameter is advantageous in a quiescent fluid but disadvantageous 
in a sufficiently strong shear flow. Locsei Il25ll previously showed that there is a 
synergy between temporal comparisons and persistence, and that when temporal 
comparisons are taken into account the drift velocity in a quiescent fluid is max- 
imised by a persistence parameter Op w 0.78, which is greater than the observed 
value of 0(p w 0.33. In that paper, it was suggested that the discrepancy between the 
observed and apparently optimal values of OCp was due to a compromise between 
transient and steady-state behaviour, with lower values of Otp favouring steady 
state behaviour. The results of the present study suggest that the observed per- 
sistence Op w 0.33 may reflect an additional trade-off between performance in a 
quiescent fluid and performance in a sheared fluid. 

The present study also shows that the elongated shape of an E. coli cell works 
to its advantage in terms of performing chemotaxis in a shear flow, particularly 
if the rotational diffusivity is zero or small. A more elongated body has a longer 
Jeffery orbit period in a shear flow, and this reduces the likelihood that the shear 
flow rotates it from facing up the gradient to down the gradient or vice-versa in 
any given run. Thus, a more elongated body maintains a positive drift velocity up 
a larger value of shear strength Q. than a less elongated body. This effect is still 
present, but less pronounced, when rotational diffusivity is accounted for, because 
the rotational diffusivity perturbs the Jeffery orbits. 

It is worth mentioning that the influence of bacterial shape and persistence of 
direction on chemotaxis in a moving fluid has been investigated in a slightly differ- 
ent context by Luchsinger et al fl6]. In their oceanographically-motivated paper 
they numerically modelled bacteria clustering around a localised nutrient source 
in a straining flow. They showed that the tendency of an elongated bacterial cell 
to align with the direction of the straining flow assists it to remain near the nutri- 
ent source. Furthermore, they showed that a 'back and forth' motion (0(p = — 1) is 
superior to a run and tumble motion (0£p = 0) in this context, because it allows the 
cell to maintain its alignment with the straining flow. Their work in conjunction 
with the present study illustrates that the 'optimal' values for chemotaxis parame- 
ters depend on context; a strategy that works well for locating a localised nutrient 
source may be different from a strategy that works well for swimming up a uni- 
form nutrient gradient. 

In conclusion, we should note that while the present work aimed to improve 
the realism of previous studies by including effects that were previously omitted, 
the model is still somewhat idealised and this should be born in mind when in- 
terpreting the results. More work, both theoretical and experimental, needs to be 
done to understand bacterial chemotaxis in realistic flow conditions. For example, 
one of the original motivations of this research was to derive the drift velocity 
and effective self-diffusivity for populations of bacteria in a general linear flow 
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(|42] |. not just a simple shear directed perpendicular to the local chemoattractant 
gradient yO. This remains a future goal. In addition, experiments indicate that in 
certain conditions boundary effects play an important role and even allow bacte- 
ria to swim upstream |21], and the fluid dynamics of this are not fully understood. 
Finally, we note that one of the strongest motivations for studying bacterial chemo- 
taxis in moving fluids is to understand the role of marine bacteria in the pelagic 
food web, and yet while it is known that marine bacteria are able to swim faster 
and react more rapidly to chemical stimuli than E. coli I27L |1], their chemotac- 
tic response has yet to be thoroughly characterised and this must be done before 
modelling can proceed on a firm footing. 

The authors are grateful for enlightening discussions with John Rallison, Ray 
Goldstein, Howard Berg, Henry Fu, Roman Stocker, Marcos, and Rachel Bearon. 
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